Orbits and reversals of a drop rolling inside a horizontal circular hydraulic jump 
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We explore the complex dynamics of a non-coalescing drop of moderate size inside a circular 
hydraulic jump of the same liquid formed on a horizontal disk. In this situation the drop is moving 
along the jump and one observes two different motions: a periodic one (it orbitates at constant 
speed) and an irregular one involving reversals of the orbital motion. Modeling the drop as a rigid 
sphere exchanging friction with liquid across a thin film of air, we recover the orbital motion and 
the internal rotation of the drop. This internal rotation is experimentally observed. 

PACS numbers: 47 55. D-, Drops and bubbles, 47.55.nb, Capillary and thermocapillary flows, 45.40.Cc, Rigid 
body and gyroscope motion, 05. 45. -a, Nonlinear dynamics and chaos. 




Usually when a drop comes into contact with the same 
liquid or a solid surface it coalesces with liquid or spreads 
on the solid. There are however exceptions, which could 
be called situations of "non wetting", when a very thin 
layer of air or vapor remains trapped between the drop 
and the substrate such as in the well known Leidenfrost 
effect pQ. In these levitation situations, the disappear- 
ance or reduction of friction with the substrate leads to 
remarkable dynamics of fluids with (nearly) no contact: 
unusual shapes (Poincare's shapes) of drops rolling down 
a plane [2], particle- wave duality of drops bouncing on 
a vibrated bath [3H5], drop motions induced by its own 
harmonic modes on vibrated viscous bath [6j [7] , chaotic 
behavior of a droplet on a soap film [3 [9], etc. 

Another case of mobile drops in a "non- wetting" situa- 
tion has been reported by Sreenivas et al. [10] and Pirat 
et al. [11 , when a drop is deposited inside a circular hy- 
draulic jump of the same liquid [T2HT5] . A thin layer of 
air is entrained underneath the drop by the supercritical 
flow of liquid feeding the jump, that prevents coalescence, 
the drop remaining trapped at the shock front with a 
strong internal rotation. In a previous letter [11] we have 
shown that this "non wetting" situation was also asso- 
ciated to remarkable dynamics in the case of a slightly 
inclined jump. In a well defined range of flow rate, a drop 
of moderate size (typically close to the capillary length) 
undergoes a gyroscopic instability, with surprising mo- 
tions along the jump perimeter, leading to oscillations 
around the lowest equilibrium position. 

In the present Letter, the authors investigate exper- 
imentally a new situation: the drop is now deposited 
inside a jump, formed on a perfectly horizontal circular 
disk, taking care to have uniform boundary conditions 
at large scale. Because of the disappearance of any refer- 
ence equilibrium position, the drop is in fact always mov- 
ing around the jump, with at least two different possible 
states: clockwise or anti clockwise regular orbital motion 
and a more complex state in which the drop "hesitates" 
between these two possible motions, leading to an irreg- 
ular behavior involving complex mechanisms of reversal. 



FIG. 1. A drop of radius a — 1 mm of silicone oil in non 
coalescence state inside a horizontal hydraulic jump of the 
same liquid (R = 5 mm). 



These reversals of motions are reminiscent of those ob- 
served in more complex hydrodynamical systems [T6HT8] , 
such as Rayleigh-Benard convection, Kolmogorov flow or 
the "dynamo" instability. They can also be observed 
in a meteorological context (quasi-biennial oscillations 
of high altitude winds) [19]. These phenomena are at- 
tracting presently a great interest from a large commu- 
nity ranging from hydrodynamics and non-linear physics 
to geophysics and meteorology. Our system constitutes 
perhaps the simplest experiment that one can build to 
observe these reversals in fluid dynamics, and in partic- 
ular without turbulence. 

A model of this new situation is proposed and pro- 
vides a reasonable description of the regular orbital mo- 
tion. As this model involves an internal rotation of the 
drop that has never been characterized, experiments were 
developed and allow us to observe and study this internal 
rotation. As predicted by our model, this rotation exists, 
but its dependence upon drop radius is more complex 
than expected from an analogy with a rigid sphere sup- 
ported by an air film under Couette flow, suggesting that 
drop shape distortions (as well a jump shape distortions) 
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FIG. 2. Spatio-temporal diagrams obtained for a drop (a = 
1.1 mm) in the periodic regime (a) (R — 4.6 mm and the 
record time is ttot — 7.4 s) and in the irregular regime (b- 
c). Case (b) (R = 5.5 mm and t to t = 28.8 s) displays a 
single reversal of the sense of rotation just above the transition 
between periodic and irregular motion, while case (c) (R — 5.9 
mm and ttot = 22.4 s) is obtained well above this one. (d) 
Phase diagram for the different kind of drop motions. (•) 
irregular motion, (♦) periodic motion and (■) drop in contact 
with the jet 



should be considered to get a better description. 

A picture of the experiment is reproduced in Fig. [T] 
A jet of silicone oil (viscosity 20 cS, surface tension 20.6 
mN/m, density 0.95) issued from a vertical tube of inter- 
nal diameter 3 mm, hits the center of a transparent glass 
disk placed 4 cm below the outlet and of radius Ro = 15 
cm. With these boundary conditions (absence of a wall 
fixing the outer height of the jump) we observed that 
the studied hydraulic jump was of type I [14 (i.e., uni- 
directional surface flow) and that the radius of the jump 
R increased with the flow rate Q while remaining very 
close to a power law R = aQ^ where a = 55 ± 10 and 
j3 = 0.77 ± 0.02 (units used for Q and R are respectively 
here m 3 .s _1 and m). A constant level tank supplied with 



liquid by a gear pump is used in order to suppress any 
flow rate pulsation. Millimeter sized drops of the same 
fluid were deposited directly inside the jump and small 
enough ones remained trapped at the shock front [TP] . 

When a drop does not touch the impinging jet (i.e. 
when the distance between the jump and the jet is larger 
than the drop width) one observes two different drop dy- 
namics: periodic and aperiodic. To characterize these 
phenomena, movies of the drop were recorded from be- 
low, through the glass plate. From the obtained frames 
the gray level evolution is extracted on a circle in order 
to obtain spatio-temporal diagrams giving the law 0(t) 
where 6 is the angular position of the drop. Three exam- 
ples are given in Fig. [2] 

(i) Fig. [2] (a) illustrates the periodic motion: the drop 
orbit ates periodically along the jump with a frequency 
ranging between a few tenths of Hz and a few Hz. This 
frequency decreases when the jump radius increases. 

(ii) Fig. [2](c) illustrates the "irregular" motion: unlike 
the periodic one this regime is not characterized by a sin- 
gle orbital period. Speed variations are observed but also 
reversals of the sense of rotation. 

(iii) In the irregular case, just beyond the transition be- 
tween periodic and irregular motion, one observes sudden 
reversals of the drop rotation sense, separating sequences 
of quasi-regular motion in opposite directions (see Fig. [2] 

GO). 

In principle, when the geometry of the experiment is 
fixed (nozzle radius, impact distance,...) one only needs 
two parameters to characterize the different regimes : the 
drop radius a and the jump radius R. A phase diagram is 
presented in Fig. [2] (d) where three different symbols are 
used for periodic and irregular motions, and for a drop in 
contact with the jet. In this last case one also observes 
complex motions (static, regular and irregular motions) 
but not well defined and with a lower drop velocity. As 
expected the dashed line separating the contact with the 
jet and the periodic motion is given by solving the equa- 
tion R = 2a + Rj where Rj is the radius of the jet that 
we find to be almost constant and equal to approximately 
1.5 mm. The continuous line separating the periodic and 
aperiodic motions is empirical but can be fitted by the 
following law R = f3a + R c where f3 = 4.4 and R c = 0.84 
mm. 

A model of the drop motion can be built as follows 
[TT] . Entrained by the radial flow of the bath, the 
drop acquires an internal rotation with a kinetic mo- 
mentum parallel to the shock front. If a perturbation 
shifts slightly the drop, the kinetic momentum conser- 
vation leads to the appearance of a radial component of 
rotation (See Fig. |3](b)). In such a situation, an active 
torque should appear, amplifying the initial perturba- 
tion and leading to a self-sustained orbital motion. To 
recover this, one models the drop as a rigid sphere of 
radius a with two contact points A and B. The situa- 
tion is described in Fig. [3] (a). We assume that the ex- 
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FIG. 3. (a) Notations and structure of the problem ex- 
pected in a radial vertical plane containing the jump central 
axis. The drop is rotating very fast because of the shear stress 
transmitted across the air film, (b) Principle of the instability 
(top view). The drop tries to keep constant its kinetic mo- 
mentum, which develops an active radial component of the 
rotation vector that tends to amplify the drop displacement. 
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changed forces at these points are simply viscous frictions 
through a sheared air film: Fa = ??a (Sa / ^a) C$ a 
Tf x gX) and Fb = Va(S B /d B )(U B -f G -t 
where r] a is the air dynamic viscosity and the quantities 
SA,S B ,dA,d B ,Tf A ,U B designate in this order : the con- 
tact surfaces, the local air layer thickness and the surface 
speed of the flow at the points A and B. G is the center of 
mass of the drop, ^ q is the speed of the center of mass 
of the drop and fl designates the angular velocity vector 
of the drop around G in the laboratory frame. 
The kinetic momentum conservation equations lead to : 

dQ r 1 L 



TA 



dt 

dVt e 1 
dQ z 1 



Mo 



dt 



r B 



— )Sl e = -u 
r B 

: + -OJ) = 

a 



'fL. 



1,Ua 
a t a 



Ub. 
tb ' 



(1) 



These equations must be coupled with the evolution 
equation for that can be deduced from the fundamental 
principles of dynamics : 

duo 2,1 1 N 2 a ^ Q z , 

at 5 ta t b 5l ta t b 

Where uj = ^ designates the orbital speed and L = 
R — a is the radius of the orbit described by the drop. 



The two characteristic times are ta — ^ 7ra ' 



r B — jFira 3 — ^- where p\ designates the mass density 
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of liquid. There are only two stationary solutions to 
these equations. A first trivial one is: 
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We have checked numerically that this solution is unsta- 
ble, in accordance with our experimental observations. 
There is also a second solution defined by: 
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which corresponds to the periodic orbital state of our 
system. Complementary experiments showed that the 
velocity Ua and U B have almost the same magnitude 
(about 0.3 m.s -1 ) and are nearly independent of the flow 
rate Q (for a typical range of 5 — 60 cm 3 .s _1 ) . It seems 
also reasonable to suppose that dA and d Bl the local air 
layer thickness at the points A and B, are independent of 
the flow rate too (they could only be fixed by the surface 
speed and by the drop properties). So one can assume 
that ta and t b are independent of the flow rate. We thus 
obtain that the orbital speed should scale as uj oc L -1 / 2 . 
This conjecture is tested in Fig. [4] by increasing the ra- 
dius (through a modification of the flow rate) for different 
drop radii. The right behavior is obtained for low jump 
radii, but a decrease for large values of L still needs to be 
explained. In addition, this model is unable to capture 
the aperiodic motion. 

The previous model is based on the hypothesis of an 
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FIG. 4. Orbital speed of the drop uj — dO/dt versus the orbit 
radius followed by the drop: L = R — a. A large range of 
drops radii is reported: 0.69 mm - 1.67 mm 
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internal rotation of the drop but this phenomenon re- 
mains to be checked. We therefore performed some ex- 
periments in order to observe this rotation by injecting 
coal particles inside the drop. Pictures of this rotation 
are reproduced in Fig. [5] The internal rotation is quite 
difficult to observe because of the drop motion inside the 
jump. To overcome this difficulty, we slightly inclined the 
jump plane (typically a few tenths of degrees), in order to 
maintain the drop globally motionless in the jump. Do- 
ing this, the solution (3) of equations (1) is selected, this 
solution being stabilized by the plate inclination [11 . As- 
suming that the speeds are proportional (or even equal), 
which seems rather reasonable, one obtains the scaling 
law: 
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Our experimental results, reproduced on Fig. [6j high- 
lighted a more complex situation : one can distinguish 
clearly two zone extremely different by increasing the 
drop radius: 

(i) For "small" drops of radii a < 1 mm (remind that 
the capillary length l c for silicone oil is about 1 mm) , in 
the limit case of the smallest radii, we recover asymptot- 
ically the - law. 

(ii) For larger drops (a > 1 mm) one can observe a 
frequency law in Qq = Ua^, where I is about 0.4 mm. 

Moreover, it is interesting to notice that all the curves 
presented here are collapsing on a single master curve 
without any rescaling. In other words the rotation fre- 
quency does not depend on flow rate. This is consistent 
with what we report above, i.e. that, in our specific situ- 
ation (boundary conditions, flow rate range, viscosity...), 
Ua and U b are nearly independent of flow rate. 

The fact that the frequency law has a cross over around 
the capillary length suggests that the hypothesis of a rigid 
sphere on a sheared gas film is too rough and that the 
drop deformations have to be taken into account, with 
also possibly some local reversal of the air flow due to 



FIG. 6. Q,g versus drop radius observed for the stabilized 
inclined static case. Different flow rates are considered but 
all the curves are collapsing on a single master curve without 
any rescaling, which suggests that the surface velocity at the 
jump could be the same for each flow rate. 



the lift applied on the drop. Indeed one can observe 
experimentally that the drops were slightly flattened by 
gravity. This is well known in the literature the drop is 
often modeled by a truncated sphere [20] , or a truncated 
sphere sustained by a non axisymetric air pocket [2T] , 

A new kind of "non wetting" dynamics for a drop in 
non coalescence state inside a horizontal hydraulic jump 
has been considered, with in particular two kinds of or- 
bital motions: periodic and irregular. A model has been 
proposed and provides good agreement for the periodic 
orbital motion. An internal rotation of the drop has been 
highlighted and validates our theories. The frequency 
law illustrates the need for further investigations on this 
fascinating object. The structure of the hydraulic jump 
has also been studied with an original mean (introducing 
drop at the shock front) and reveals some unexpected 
observations about the local surface speed at the jump 
(does it depend or not upon flow rate ?) that would 
deserve further investigations. 
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